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A  Numerical  Analysis  of  the  Effect  of  Inhomogeneity  and  Wind 
Speed  Shear  on  the  Vertical  Dispersion  in  the  Mesoscale  Range 

Lei  Xiaoen  and  Yan  Bangliang* 

(Institute  of  Atmospheric  Physics,  Academia  Sinica) 


Abstract 

According  to  Expression  (10)  in  Reference  [1] ,  the  effect 
of  inhomogeneity  and  wind  speed  shear  on  the  vertical  dispersion 
in  the  mesoscale  range  (10-lQ0km)  was  analyzed  numerically. 

Some  beneficial  results  were  obtained. 

In  this  paper,  starting  from  the  correlation  between  the 
vertical  dispersion  parameter  a  and  x  derived  in  [1]  (the 

Z 

meanings  and  parametric  values  of  all  the  symbols  used  in  this 
paper  are  identical  to  those  in  [1] ,  and  in  general  will  not  be 
explained  again) ,  the  effect  of  the  variation  of  the  effective 
roughness  length  zQ  with  distance  (i.e. ,  inhomogeneity)  and 
various  wind  profiles  on  a  was  analyzed  numerically.  It  pro- 
vided  a  method  to  further  resolve  practical  problems  (such  as 
the  evaluation  of  the  atmospheric  quality  in  the  Beijing-Tianjin- 
Bohai  area) . 

I.  Effect  of  Wind  Shear  Speed  on  a 

z 

From  the  equation  (10)  in  [1]  one  can  see  that  the  effect 

of  the  variation  of  the  wind  speed  with  the  altitude  on  o  is  an 

z 

important  factor.  In  the  following,  besides  using  equation  (18) 
(profile  II)  given  in  [1] ,  we  also  chose  two  other  types  of  wind 
profile  to  carry  out  a  comparative  numerical  analysis. 

In  order  to  describe  the  upper  half  of  the  PBL  better,  a 
two  layer  model  such  as  the  one  used  in  Reference  L2]to  analyze  the 
Eikman  laminar  wind  speed  profile  (i.e.,  considering  a  transi¬ 
tion  layer  z^,  u*  is  a  constant  below  z^,and  KQ  is  a  constant 
above  z^)  was  adopted.  Above  z^ ,  all  the  way  to  an  altitude  h, 


— (l+,')V'?/i*'V(7— r,» 


(i) 


where 


«*«■«+*» 


(2) 


Based  on  (1)  and  (2)  ,  a  complex  number  calculation  was  conducted. 
Finally,  we  obtained 


u—  (?  cos  o, + cos  -(*-»,)/  J 


(3) 


Above  the  constant  stress  layer,  equation  (3)  (Profile  III)  was 
more  adequate  than  II.  Below  z^,  it  was  possible  to  apply  the 
near  ground  layer  law.  When  z^  is  very  low,  equation  (3)  can 
approximately  describe  the  wind  speed  distribution  characteris¬ 
tics  of  the  near  ground  layer. 


♦Comrade  Yaun  Suzhen  participated  in  part  of  the  calculation  and 
analysis  work. 

The  power  series  wind  profile  has  been  widely  used  by  vari-  /75 
ous  industrial  departments.  For  comparison,  the  following. wind 
profile ^  was  also  chosen  (Profile  I)  in  the  analysis:  (4) 

During  the  computation,  r  was  chosen  to  be  0.05  in  unstable  con- 
[41 

ditions,  r  =  0.1  in  neutral  conditions,  and  r  =  0.4  in  stable 
conditions . 

Which  one  of  the  wind  profiles  I,  II,  and  III  agrees  more 

with  the  reality?  The  dimensionless  wind  profiles  (u/u2Q)  with 

various  stabilities  were  plotted  in  Figure  1.  For  comparison 

purposes,  some  of  the  wind  profiles  observed  on  a  430m  tower  in 

[41 

the  United  States  and  a  320m  tower  in  Beijing  |_5]  were  also  plot¬ 
ted  in  dots  in  the  figure. 

From  Figure  1  one  can  see  that  the  distribution  of  I  is 
different  from  II  and  III .  II  and  III  have  the  same  distribu¬ 
tion  shape.  '  Among  the  three,  II  has  the  largest  vertical  wind 
shear.  The  wind  shear  of  I  is  close  to  that  of  III.  A  com¬ 
parison  with  the  wind  actually  measured  shows  that  III  is  more 
close  tc  reality. 

The  expressions  of  I,  II,  and  III  were  substituted  into  the 

relation  between  oz  and  x.  Numerical  computations  were  carried 

out  on  a  TQ-6  machine.  The  u*  values  are  given  in  Table  1. 

Finally,  the  magnitudes  of  a  were  obtained  under  different  wind 

profile  conditions.  In  order  to  examine  the  difference  among 

them,  the  ratios  of  o.  obtained  under  the  three  conditions  are 

z 

given  in  Table  2 . 


r  0.355 
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0.32  i  0.28 


0.244 


0.16 


Table  1.  Values  of  u*  at  various  stabilities^”9^ 

1.  type  of  stability 
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Table  2.  Ratios  of  az  derived  with  various  wind  profiles* 

1 .  type  B 

2 .  type  D 

3 .  type  E 

4.  *in  which  a  (I) ,  o  (II)  ,  and  o  (III)  represent  the 

results  obtained  using  Profiles  I,  11,  and  III  in 
the  computation,  respectively. 
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The  ratio  of  a  and  the  profile  comparative  analysis  in 
z 

Figure  1  showed  that: 

1.  In  a  neutral  condition,  the  wind  shears  of  II  and  III 
can  differ  by  as  large  as  nearly  20  times.  The  a  ratio  thus 
produced  has  a  maximum  of  1.39  (in  the  10-100  km  range).  Its 
average  is  1.34  i  0.082.  The  shapes  of  Profiles  II  and  I 
are  not  the  same.  The  wind  shear  difference  of  the  two  differs 
by  as  much  as  19  times . 


•  A»  »• 


Figure  1.  The  Average  Wind  Speed  Shear  Profile 

1.  B.C.D.E:  stability  type 

2.  I,  II,  III:  profile  type 

3.  observation  results  obtained  at  a  320m  tower 

4.  dotted  line:  observation  results  obtained  at  a  430m  tower  in  the 
United  States 

3.  altitude 


/■.  r . 


The  maximum  az  ratio  is  1.45.  Its  average  is  1.35  ±  0.099.  /77 

Comparing  III  to  I ,  despite  the  fact  that  the  profiles  type  are 

not  the  same,  the  wind  speed  shear  difference  is  very  small. 

The  difference  in  a  is  also  small.  The  maximum  ration  is  only  1.07. 

2.  From  the  variation  of  a  derived  from  both  II  and  III 

z 

with  stability,  the  difference  of  a  becomes  greater  and  greater 
as  the  type  changes  from  A  to  D .  This  is  primarily  because  the 
wind  shear  of  II  is  increased  gradually  from  type  A  to  D.  There¬ 
fore,  a  gradually  becomes  smaller.  The  wind  shear  of  III,  how- 
ever,  does  not  vary  significantly  with  stability.  The  same 
trend  of  variation  also  exists  from  type  F  to  D.  However,  it  is 
not  as  obvious  as  under  an  unstable  situation.  The  pattern  of 
variation  of  the  effect  of  both  II  and  I  on  az  with  stability  is 
the  same  as  that  of  II  and  III.  Comparing  III  to  I,  due  to  the 
fact  that  the  variation  of  the  wind  shear  with  stability  is  not 
very  obvious  for  both  cases ,  the  wind  shear  is  very  close . 

Although  the  distribution  shapes  are  different,  the  difference 
of  the  effect  of  the  two  on  az  is  always  within  10%. 

II.  Effect  of  Inhomogeneity  on  az 

Due  to  the  effect  of  roughness  on  dispersion,  the  discus¬ 
sion  in  Reference  [1]  used  an  effective  roughness  length  to  ex¬ 
press  the  average  value  of  a  homogeneous  geomorphology  and  an 
inhomogeneous  geomorphology  after  an  abrupt  change  takes  place. 

This  is  feasible  in  the  estimation  of  the  average  condition. 

However,  in  reality  the  roughness  factor  is  frequently  inhomo¬ 
geneous.  It  is  a  problem  of  our  great  concern  to  know  how  large 
the  effect  of  inhomogeneity  is.  For  simplicity,  but  not  losing 
its  generality,  Zq  is  expressed  as 

■  «,(i)-cos(JILx')+l.QQ05  (5) 

>  \  / 

For  convenience,  — ! -  was  chosen  to  be  1 .  The  value  of  zn  given 

u  -  u 

by  equation  (5)  in  fact  is  varying  periodically  with  x  between 

0 .0005-2 .0005m.  With  regard  to  such  a  roughness  distribution, 

the  effective  roughness  length  is  approximately  lm. 
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Figure  2.  Effect  of  zQ  (X)  on  a 

1.  the  solid  line  is  for  zQ=0.0  1  m 

2.  the  dotted  line  is  for  z»=(x) 

3.  Except  for  curve  1,  which  is  for  z0=5m,  other 
curves  are  for  the  situation  zQ=l. 

In  order  to  compare  the  calculated  results  to  the  value  of 
<j2  when  Zg  is  a  constant,  it  is  assumed  that  the  inhomogeneity 
in  the  form  of  equation  (5)  begins  at  3km  away  from  the  source 
in  the  calculation.  Equation  (5)  is  substituted  into  the  cor¬ 
relation  between  a  and  x  to  carry  out  the  numerical  calcula¬ 
tion.  The  results,  together  with  the  values  obtained  when  z^  - 
is  a  constant,  are  given  in  Figure  2.  In  the  meantime,  the 

mean,  maximum,  and  minimum  values  of  a  (x)/a  (in  the  range  of 

z0  zo 

5-100km)  in  various  types  of  stability  are  tabulated  in  Table  3. 


O  » 


r 


* 


y 


> 


I'*,  (i/'ul 

<»V'| 


,  **,  (£i/**si 

2  <'JS 


(A/*l 


l.68±0.361 

Q95±0.023 

0.55+0.10T 


UI  +  0.Q23 
0.08±  0.008 
0.88  +  0  028 


1.59x0108 
159  +  0.0T3 
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Table  3.  The  Values  of  a_  (x)/a„  * 

z0  z0 

1 .  Type  A 

2 .  Type  D 

3 .  Type  F 

4 .  maximum 

5 .  mean 

6 .  minimum 

7.  *The  subscript  of  a  represents  the  roughness  length 

(m) .  zq(x)  represents  inhomogeneous  situations 
and  zQ  represents  homogeneous  situations. 

From  Table  3  and  Figure  2  one  can  see  that: 

1 .  When  the  pollutant  abruptly  enters  an  inhomogeneous 
terain  in  the  form  of  equation  (5)  after  traveling  originally  in 
a  flat  and  homogeneous  terrain,  the  a2  versus  x  curve  shows  an 
apparent  change  in  the  corresponding  region  to  form  a  turning 
point.  Due  to  the  fact  that  the  average  roughness  of  the  in¬ 
homogeneous  geomorphology  is  greater  than  0.01m,  the  value  of  a 

z 

is  larger  than  that  corresponding  to  a  flat  homogeneous  terrain. 


The  mean  value  of  a. 


(x)/a0  01  becomes  larger  as  it  varies  from 


D  to  F.  From  Type  D^to  A,  the  variation  trend  is  not  obvious. 

The  average  ratio  of  Type  D  is  the  smallest. 

2.  When  the  pollutant  abruptly  enters  an  inhomogeneous 

region  in  the  form  of  equation  (5)  from  a  homogeneous  terrain 

whose  roughness  is  lm,  the  results  indicate  that  o  ,-,/a.sl 

z  ix;  i 

from  Type  D  to  A  in  neutral  situations.  It  is  basically  very 
close  to  the  homogeneous  situation  in  which  zn  =  lm.  However, 


r*  rvr 


the  stable  situation  is  not  quite  the  same,  a  is  appar- 

Zq  (X  ) 

ently  larger  than  when  Zg=lm.  This  phenomenon  cells  us  that, 
with  regard  to  unstable  and  neutral  conditions,  it  is  as  ef¬ 
fective  to  use  a  homogeneous  terrain  with  an  effective  rough¬ 
ness  Zq  to  represent  an  inhomogeneous  rough  terrain  with 
average  roughness  zQ  to  investigate  the  effect  on  c  .  The 
stable  situations  are  quite  different.  This  inhomogeneity 

could  also  cause  a  turbulent  disturbance ,  which  caused  a  to  be 

z 

larger  than  the  average  roughness  length  situation.  This  fur¬ 
ther  explains  the  importance  of  the  inhomogeneity  of  the  terrain 
with  respect  to  stability. 

3 .  When  the  pollutant  abruptly  enters  an  inhomogeneous 
region  in  the  form  of  equation  (5)  from  a  terrain  whose  rough¬ 
ness  is  greater  than  lm,  the  results  indicate  that  inhomogeneity 
would  cause  a  to  produce  a  turning  point  under  unstable  and 

neutral  conditions.  But,  the  value  of  a  is  smaller  than  that 

z 

obtained  when  zQ  is  5m.  This  difference  decreases  with  in¬ 
creasing  distance.  It  is  apparently  different,  however,  under 
stable  conditions.  Regardless  of  the  initial  homogeneous  top¬ 
ology,  the  value  of  a  is  always  larger  than  that  obtained  when 

z 

zQ  does  not  vary  due  to  inhomogeneity.  This  demonstrated  the 
inhomogeneity  of  the  roughness  is  purely  mechanical  in  neutral 
and  unstable  situations .  As  for  stable  situations ,  in  addition 
to  a  mechanical  action,  it  indirectly  intensifies  the  turbulence 
This  mechanism  is  yet  to  be  further  understood. 

Ill .  Conclusions 

Through  the  above  analysis,  one  can  obtain  that: 

1 .  In  the  mesoscale  range ,  the  distribution  of  wind  speed 

with  altitude  and  the  vertical  wind  speed  shear  have  a  greater 

effect  on  a  .  The  magnitude  of  the  wind  shear  is  more  important 
z 

than  the  profile  type. 

2.  Inhomogeneity  could  cause  a  to  form  a  turning  point 

M 

at  the  place  where  the  terrain  changes  abruptly. 

3 .  It  is  possible  to  use  an  effective  roughness  length  to 
replace  the  inhomogeneity  in  neutral  and  unstable  situations. 


However,  in  stable  situations,  it  would  cause  a  large  error. 

Received  on  January  15,  1981. 
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